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Abstract 

We present a detailed, all-order study of gluon mass generation within the massless bound- 



o: 

state formalism, which constitutes the general framework for the systematic implementation of the 
Schwinger mechanism in non-Abelian gauge theories. The main ingredient of this formalism is 

. 

the dynamical formation of bound-states with vanishing mass, which give rise to effective vertices 

(N" 

containing massless poles; these latter vertices, in turn, trigger the Schwinger mechanism, and 
allow for the gauge- invariant generation of an effective gluon mass. This particular approach has 
the conceptual advantage of relating the gluon mass directly to quantities that are intrinsic to the 



bound-state formation itself, such as the "transition amplitude" and the corresponding "bound- 
state wave- function" . As a result, the dynamical evolution of the gluon mass is largely determined 



by a Bethe-Salpeter equation that controls the dynamics of the relevant wave- function, rather than 



the Schwinger-Dyson equation of the gluon propagator, as happens in the standard treatment. 
The precise structure and field-theoretic properties of the transition amplitude are scrutinized in 
a variety of independent ways. In particular, a parallel study within the linear covariant (Landau) 
gauge and the background-field method reveals that a powerful identity, known to be valid at 
the level of conventional Green's functions, relates also the background and quantum transition 
amplitudes. Despite the differences in the ingredients and terminology employed, the massless 
bound-state formalism is absolutely equivalent to the standard approach based on Schwinger-Dyson 
equations. In fact, a set of powerful relations allow one to demonstrate the exact coincidence of 
the integral equations governing the momentum evolution of the gluon mass in both frameworks. 

PACS numbers: 12.38.Aw, 12.38.Lg, 14.70.Dj 
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I. INTRODUCTION 



The dynamical generation of an effective (momentum-dependent) gluon mass in pure 
Yang-Mills theories [l-3] (and eventually in QCD) has attracted considerable attention in 
recent years, because it furnishes a cogent theoretical explanation for a large number of im- 
portant lattice findings j4-ll|. This inherently nonperturbative phenomenon is usually stud- 
ied within the formal machinery of the Schwinger-Dyson equations (SDEs), supplemented by 
a set of fundamental guiding principles, which enable the emergence of "massive" solutions, 



while preserving intact the gauge invariance of the theory 
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The main theoretical concept underlying the existence of such solutions is the Schwinger 



mechanism 
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151 ] . The basic observation is that if the "vacuum polarization" acquires a 



pole at zero momentum transfer, then the vector meson becomes massive, even if the gauge 
symmetry forbids a mass at the level of the fundamental Lagrangian, as happens in the 
case of gauge theories. In the absence of elementary scalar fields, the origin of the afore- 
mentioned poles must be purely nonperturbative: for sufficiently strong binding, the mass 
of certain (colored) bound states may be reduced to zero 16|-[20|]- In addition to triggering 
the Schwinger mechanism, these bound-state poles act as composite, longitudinally coupled 
Nambu-Goldstone bosons, enforcing the gauge invariance of the theory in the presence of 
gauge boson masses. Every such Goldstone-like scalar, "absorbed" by a gluon in order to 
acquire a mass, is expected to actually cancel out of the S'-matrix against other massless 
poles or due to current conservation. 

The concrete realization of the above general scenario in Yang-Mills theories proceeds by 
appealing to the existence of a special type of nonperturbative vertices, to be generically 
denoted by V (also referred to as "pole vertices"), which essentially transmit the effects of 
the composite excitations to the dynamical equations governing the various off-shell Green's 
functions. Specifically, the inclusion of these vertices (i) enables the SDE of the gluon 
propagator to admit massive solutions, and (m) guarantees that the Ward identities (Wis) 
and the Slavnov- Taylor identities (STIs) of the theory maintain exactly the same form before 
and after gluon mass generation; a particular consequence of this important property is the 
exact transversality of the resulting (massive) gluon self-energy. 

As has been demonstrated in a recent work [21], the full dynamical equation of the 
gluon mass in the Landau gauge may be derived from the SDE of the gluon propagator by 
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postulating the existence of the pole vertices, and employing the Wis and STIs they satisfy, 
together with their totally longitudinal nature (for related studies in the Coulomb gauge, 



sec 
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25] ). The resulting integral equation makes no reference to the closed form of these 



vertices, nor to the actual dynamical mechanism responsible for their existence. However, 
the details of the actual dynamical formation of the bound-state poles, and subsequently 
of the pole vertices, as well as their exact closed form, transcend the concrete purpose of 
obtaining a particular version of the gluon mass equation, and are of paramount importance 
for the self-consistency of the entire gluon mass generation scenario. In fact, it would be 
highly desirable to establish a precise quantitative connection between the fundamental 
ingredients composing these vertices and the gluon mass itself. 

The purpose of the present work is to dissect the pole vertices and scrutinize the field- 
theoretic properties of their constituents, within the context of the "massless bound-state 
formalism", first introduced in some early seminal contributions to this subject 
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20|, 



and further developed in 26]. The final outcome of this analysis is an alternative, but 
completely equivalent, description of the dynamical gluon mass in terms of quantities ap- 
pearing naturally in the physics of bound states, such as the "transition amplitude" and 
the "bound-state wave function" . This new description makes manifest some of the salient 
physical properties of the gluon mass (e.g., positive-definiteness), and provides, in addition, 
a decisive confirmation of the self-consistency of the concepts and methodology employed. 

Let us next present the general outline of the article, introducing some of the basic 
concepts, and commenting on the logical connections and delicate interplay between the 
various sections. 

The starting point of our considerations is a brief review of the importance of the pole 
vertices for obtaining massive solutions out of the SDE for the gluon propagator, in a gauge 
invariant way, i.e., preserving the form of the fundamental STIs of the theory (Section ITT|) . 

Within the massless bound-state formalism, the pole vertices are composed of three fun- 
damental ingredients, (i) The nonperturbative transition amplitude, to be denoted by I(q 2 ), 
which connects a single gluon to the massless excitation, (ii) The scalar massless excitation, 
whose propagator furnishes the pole i/q 2 , and (Hi) a set of "proper vertex functions" [l7| (or 
"bound-state wave functions"), to be generically denoted by B, (with appropriate Lorentz 
and color indices), which connect the massless excitation to a number of gluons and/or 
ghosts. The quantity I(q 2 ) is universal, in the sense that it appears in all possible pole ver- 
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tices. Furthermore, it admits its own diagrammatic representation, which, in turn, involves 
the functions B. As a result, the dependence of the pole vertices on the Bs is quadratic 
(Section Ell). 

When inserted into the SDE for the gluon propagator, and all diagrams are kept (no 
truncation), the special structure of the pole vertices allows one to obtain a very concise re- 
lation between the gluon mass and the square of the transition amplitude, given in Eq. (14. 5p . 
This relation demonstrates that, unless I(q 2 ) vanishes identically, the gluon mass obtained 
is positive-definite (Section HV|) . 

If the above construction is repeated within the combined framework of the pinch tech- 
nique (PT) 



scheme' : 
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3l| and the background field method (BFM) |32j (known as the "PT-BFM 



3-35]), the relevant quantity to consider is the transition amplitude between a 
background gluon and the massless excitation, to be denoted by I{q 2 ) (all other ingredients 
remain identical). The mass of the propagator connecting a background (B) and a quantum 
gluon (Q) can then be expressed as the product of I(q 2 ) and I(q 2 ); an equivalent relation may 
also be obtained from the background-background propagator (BB). The additional fact 
that the PT-BFM gluon propagators (QB and BB) are related to the conventional one (QQ) 



by a set of powerful identities (known as Background-Quantum identities (BQIs) 36|, |37J), 
allows one finally to relate I(q 2 ) and I{q 2 ) as shown in Eq. fl 5 . 1 2 [) (Section |V|) . 

Evidently, Eq. (" 15 . 1 2[) has emerged as a self-consistency requirement between two different 
formulations of the SDEs (conventional and PT-BFM), which, in their untruncated version, 
must furnish the same physics. It would be very important, however, to establish the validity 
of Eq. (I5.12p in a more direct, explicit way, by operating at the level of the pole vertices 
themselves, where I{q 2 ) and I{q 2 ) make their primary appearance. To that end, we will 
carry out the explicit construction of the three-gluon pole vertex, both in the conventional 
and the BFM formalism, using as a sole input the WI and/or STIs they satisfy, and their 
totally longitudinal nature (subsection IVI Al) . The comparison of the two, after judicious 
identification of the parts that contribute to the gluon mass equation, reproduces precisely 
Eq. (jgZEZD (subsection MM- 

It turns out that an even more fundamental derivation of Eq. (15. 1 2[) may be devised, which 
takes one back to the underpinnings of the PT-BFM connection: the BQIs, which are for- 



mally obtained within the Batalin-Vilkovisky (BV) formalism 
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can be alternatively 



derived through the diagrammatic rearrangements implemented by the PT (Section IVIII) . In 
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the case of the SDE series containing regular (fully dressed) vertices, such a diagrammatic 
derivation amounts finally to the demonstration of the nonperturbative PT-BFM equiva- 
lence 35] . From the operational point of view, the standard PT construction boils down to 



the judicious exploitation of the rearrangements produced when certain longitudinal (pinch- 
ing) momenta trigger the STIs satisfied by the aforementioned vertices; the required STIs 
are known from the formal machinery of the BV, or alternative formalisms. A priori, the 
implementation of this procedure at the level of the diagrammatic representation of the I(q 2 ) 
is thwarted by the fact that the pinching momenta will act on the B vertices (introduced 
above), whose STIs, however, are not formally known. Nonetheless, the explicit construc- 
tion of the pole vertex, and the subsequent line of reasoning, furnish precisely the missing 
STI, and make the PT-driven diagrammatic construction possible (subsection IVH Ap . An 
interesting by-product of this construction is the derivation of an integral constraint between 
the S-vertex functions containing ghost legs. When combined with the recent lattice find- 
ings on the infrared behavior of the ghost propagator, this constraint strongly suggests the 
individual vanishing of all such ghost vertex functions (subsection IVH 61 . 

The diagrammatic evaluation of the transition amplitude furnishes Eq. (I8.15p . which 
expresses this fundamental ingredient of the bound-state formalism in terms of a double 
integral containing solely the vertex B^ v . Even though this vertex function is not known 
for general momenta, the powerful relation of Eq. ( I8.18p . first derived in [26], relates its 
g^v form factor to the gluon mass. This relation, in turn, allows one to recover the mass 
equation derived in [21], following a completely different methodology and formalism. The 
resulting exact coincidence reveals an impressive complementarity between two formally 
distinct methods (subsection IVIH A~|) . 

The equivalence of the two formalisms established in the previous subsection is exact only 
when no approximations are employed in either one. However, in practice approximations 
must be carried out, and as a result, in general, due to the difference of the intermediate 
steps, a modified momentum dependence for the gluon masses will be obtained. In subsection 
IVIIIBI we outline the general computational procedure that must be followed within the 
massless bound-state formalism in order to determine the dynamical gluon mass. The nature 
of this procedure is a-priori dissimilar to that employed in 2l|; thus, whereas an SDE 
equation is solved in the latter case, now the main dynamical ingredient is the Bethe-Salpeter 
equation (BSE) that controls the evolution of the relevant bound-state wave function. In 
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the qualitative discussion presented here, we focus on the difference in the momentum- 
dependence of the gluon mass that one observes at the lowest level of approximation within 
both formalisms. 

Finally, our conclusions and discussion are presented in Section IIXI . 



II. GETTING MASSIVE SOLUTIONS FROM THE GLUON SDE 

In this section we review the general principles that allow the generation of massive 
solutions out of the gluon SDE, and study in detail the general structure of the special pole 



vertices that trigger this effect (for different approaches, see, e.g., |40|-|55|). 



A. General principles 

The full gluon propagator A^(g) = 5 ab /S.^ u (q) in the Landau gauge is defined as 



A^(g) = -^P^(g)A(g 2 ), (2.1) 



where 



PM = 9r- S f L , (2-2) 

is the usual transverse projector, and the scalar cofactor A(g 2 ) is related to the (all-order) 
gluon self-energy U^q) = P^(q)U(q 2 ) through 

A-V)=g 2 + ffl(g 2 ). (2.3) 
Alternatively, one may introduce the inverse of the gluon dressing function, J(q 2 ), defined 



as 



56] 



A-V) = g 2 J(g 2 ). (2.4) 

Let us consider the SDE for the gluon propagator in the conventional formulation of Yang- 
Mills (i.e., linear gauges) 

&-\q 2 )P^{q) = q 2 P„(q) + Hl^q), (2.5) 
with the self-energy given by 

5 

n^(g) = X)( a i)^ 5 (2-6) 

i=l 
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FIG. 1: SDE satisfied by the conventional gluon self-energy, with two quantum gluons (Q) entering. 
White (black) circles denote fully-dressed propagators (vertices). 

where the diagrams (a*) are shown in Fig. [TJ Note that the full (untruncated) self-energy is 
transverse, namely 

5 

g"lW<z) = Xy«U = 0- (2-7) 



i=l 



It turns out that, with the Schwinger mechanism turned off (i. e., in the absence of massless 
poles in the vertices) the SDE leads to the conclusion that A _1 (0) = 0, namely the absence 
of massive solutions. 

In order to obtain massive solutions out of the above SDE, and preserve, at the same 
time, the gauge invariance intact, one must carry out the crucial substitution 

r _ ► r = r m + v, (2.8) 

to all fully-dressed interaction vertices appearing in Eq. ( 12. 5ft . The main characteristic of 
the vertices V, which sharply differentiates them from ordinary vertex contributions, is that 
they contain massless poles, originating from the contributions of bound-state excitations. 
Such dynamically generated poles are to be clearly distinguished from poles related to ordi- 
nary massless propagators, associated with elementary fields in the original Lagrangian. In 
addition, they are completely longitudinally coupled, i.e., they satisfy conditions of the type 
(for the case of the three-gluon vertex) 

P a ' a (q)Pv'^r)P^(p)V a/ ,, u/ (q,r,p) = 0. (2.9) 

As for the vertices T m , they are given by the same graphs as the T before, but with gluon prop- 
agators replaced by massive ones [see Eq. (I2.12p ], implementing simultaneously Eq. f |2.8j) . 



The new (massive) self-energy is then given by 
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IW<Z) = J>0m- (2-10) 



where the "prime" indicates that the various fully-dressed vertices appearing inside the corre- 
sponding diagrams of the gluon self-energy have been replaced by their primed counterparts, 
as in Eq. (|2.8)) . It is important to emphasize that, since the above replacement maintains 
the STIs of the theory unaltered, the transversality of the massive self-energy persists, i.e., 

5 
i=l 

The appearance of massive solutions amounts effectively to the change (in Minkowski 
space) 

A" V) = q 2 J(q 2 ) —> A" V) = q 2 J m {q 2 ) ~ m 2 (q 2 ), (2.12) 

with m 2 (0) 7^ (of course, in Euclidean space, one must find A~ 1 (0) > 0). The subscript 
"m" in J m indicates that effectively one has now a mass inside the various expressions: for 
example, whereas perturbatively J{q 2 ) ~ lng 2 , after dynamical gluon mass generation has 
taken place, one has J m (q 2 ) ~ ln(g 2 + m 2 ). 

The actual evaluation of the relevant diagrams may be carried out by appealing to the 
basic global features of the V vertices, as deduced from the STIs and their complete longi- 
tudinality 21] . The final upshot is that the SDE may be schematically cast into the form 
(Minkowski space) 

q 2 J m {q 2 ) ~ m 2 {q 2 ) = q 2 [l + K^q 2 , m 2 , A m )] + K 2 {q 2 , m 2 , A m ), (2.13) 

such that g 2 /Ci(g 2 , m 2 , A m ) — > 0, as q 2 — )> 0, whereas JC 2 (q 2 , m 2 , A m ) ^ in the same limit, 
precisely because it includes the term 1/g 2 contained inside V afiu (q,r,p). This form, in turn, 
gives rise to two coupled integral equations, an inhomogeneous equation for J m (q 2 ), and a 
homogeneous one for m 2 (q 2 ) (the latter is usually referred to as the "mass equation"), of 
the generic type 

J m (q 2 ) = 1+ / /d(g 2 ,m 2 ,A m ), (2.14) 

m 2 (q 2 ) = - [ /C 2 (g 2 ,m 2 ,A m ). (2.15) 
Jk 

Physically meaningful (i.e., positive definite and monotonically decreasing) solutions to an 



approximate version of the full mass equation have been recently presented in 



2l|. 



B. Structure of the vertices V 



We will next consider the decomposition of the vertices V that emerges in a natural way, 
if one employs as a criterion the effect that the various components of V may have on the 
gluon SDE. For concreteness we will focus on the case V afJiU ; however, all basic arguments 
may be straightforwardly extended to any vertex of the type V. 

Since the main function of the vertex V atiu is to generate a mass term when inserted into 
the graph (ai) of Fig. [Q, it must contain components that do not vanish as q — > 0. To study 
this point in more detail, we first separate V into two distinct parts, namely 

Vafiu{q,r,p) = U aflv (q,r,p) + R aiMU (q,r,p) , (2.16) 

defined as follows. U is the part of V that has its Lorentz index a saturated by the momen- 
tum q; thus, it contains necessarily the explicit g-channel massless excitation, namely the 
1/q 2 poles. It assumes the general form 

U aiiV {q,r,p) = ^C^,(g,r,p), (2.17) 

where, due to Bose symmetry under the exchange r p, \i «->• u, C^ v {q, r, p) must satisfy 
C„^(q,p,r) = -C^(q,r,p); as a result, C^(0, -p,p) = 0. 

The term R contains everything else; in particular, the massless excitations in the other 
two kinematic channels, namely 1/r 2 and 1/p 2 (but not 1/q 2 ) are assigned to R. Thus, for 
example, terms of the form q a 9ij,u and qat^Pv are assigned to U alJiU , while terms of the type 
PvQcui and r^g aL , belong to R. Evidently, P^'^(r)P u ' u (p)R a ^ I/ /(q,r,p) = 0. In addition, again 
due to Bose symmetry, we have that R afJ ,v(0, ~P,p) = 0. 

Consider now the individual effect that U and R have when inserted into the gluon SDE, 
specifically the graph (ai) of Fig. [TJ Of course, it is expected, on general grounds, that R 
should not generate mass terms, because it does not contain poles of the type 1/q 2 ] this is 
indeed what happens. If we work in the Landau gauge, any contribution from R vanishes 
for a simple kinematic reason, namely the above transversality condition satisfied by R. 
Away from the Landau gauge, R does not contribute to the mass equation Eq. (12.151) either, 
because R atJ ,v(0, —p,p) = and there is no 1/q 2 term that could compensate this. Thus, R 
seems to contributes, in a natural way, to the equation for J m (q 2 ). 

In fact, it is relatively straightforward to establish that, in the limit q — > 0, the SDE 
contribution generated by R vanishes as 0(q 2 ) (or faster). Indeed, let us set p = k, the 
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virtual integration momentum in the SDE diagram; given that R af iv{Q, —k, k) = 0, a Taylor 
expansion of R around q 2 = gives (suppressing indices) 



Now, the first term is odd in k, and therefore the corresponding integral vanishes. As 
a result, the term 2(qk)R'(k 2 ) must be multiplied by another term proportional to (qk), 
coming from the rest of the terms (propagators, vertices, etc) appearing in the integrand. 
Thus, the resulting contribution is of order 0(q 2 ) (or higher), as announced. The importance 
of this property is related to the fact that, in order to arrive at Eq. (I2.14p . one must pull 
out of the corresponding integral equation a factor of q 2 \ had the integrand vanished slower 
than 0(q 2 ), one would get a divergent contribution for J m (0), which would be physically 
unacceptable, given that J m (0) is inversely proportional to the infrared finite QCD effective 
charge [see, e.g., (5^], and references therein]. 

Therefore, the only term that can contribute to the gluon mass equation is U a ^ u (q,r,p)] 
the precise contribution will depend, of course, on the exact form and behaviour of the 
cofactor C^(g, r, p), as q — > 0. It is clear for instance, that if C^ u (q, r, p) contains terms that 
behave as 0(q 1+c ), with c > 0, as q — >■ 0, then these terms could not possibly trigger the 
Schwinger mechanism, because the effect of the pole would be counteracted by the positive 
powers of q. On the other hand, terms that vanish as C(g 1_c ), would give rise to divergent 
results; however, this latter possibility does not occur, again due to the Bose symmetry of 
V with respect to p <->■ r. Thus, the only terms of C^ u relevant for gluon mass generation 
are those that vanish as 0(q) (c = 0). 

These observations motivate the separation of U atJ-u (q,r,p) into two parts, one that be- 
haves as a constant, O(q ), thus contributing to the mass equation (to be denoted by U a ^ u ), 
and one that vanishes as 0(q c ), c > 0, and can be naturally reassigned to R (to be denoted 
by W' ). In fact, due to the same reason explained above, namely, the absence of divergent 
contributions at the level of Eq. (I2.14p . we must have that c > 1; actually, the explicit 



construction presented in Section |Vl] [e.g. Eq. (I6.24p and ensuing discussion], reveals that 
c = 1. Thus, 



R{q, -q -k,k) = 2{qk)R'{k 2 ) + 0{q 2 ) . 



(2.18) 




(2.19) 
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Then, setting 

fc afi u(q,r,p) = R aflv (q,r,p) + U^ u {q,r,p) (2.20) 
we arrive at the final separation 

V a ^(q,r,p) = U ai i V (q,r,p) + lZ af , u (q,r,p^ (2.21) 

m 2 (g 2 ) J m (g2) 

where the curly brackets below each term indicate to which equation [Eq. f !2.14p or Eq. (I2.15P 
each term will contribute. 

III. MASSLESS BOUND-STATE FORMALISM 

Whereas in the SDE approach outlined in the previous section one relies predominantly 
on the global properties of the vertices V, within the massless bound-state formalism one 
takes, instead, a closer look at the field-theoretic composition of these vertices, establishing 
fundamental relations between their internal ingredients and the gluon mass. This becomes 
possible thanks to the key observation that, since the fully dressed vertices appearing in 
the diagrams of Fig. [TJ are themselves governed by their own SDEs, the appearance of such 
massless poles must be associated with very concrete modifications in the various structures 
composing them. 

Let us begin by recalling that, in general, when setting up the usual SDE for any fully- 
dressed vertex contained in Eq. (12.51) . a particular field (leg) is singled out, and is connected 
to the various multiparticle kernels through all elementary vertices of the theory involving 
this field (leg). The remaining legs enter into the various diagrams through the aforemen- 




FIG. 2: The SDE for the Q 3 vertex T aiMU (q,r,p). Gray blobs denote the conventional 1PI (with 
respect to vertical cuts) multiparticle kernels. 
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(A) 




(a) (6) 




FIG. 3: (A) A diagram not included into the standard kernel. (B) The kernel with the Schwinger 
mechanism turned on: in addition to the "regular part" (b'i) (gray striated), the massless excitation 
in the (/-channel (b' 2 ) is added. (C) The part (b'i) is obtained from the original (gray) kernel (6) by 
inserting massive gluon propagators into its diagrams, and carrying out the substitution Eq. fj2.8|) 
in the fully-dressed vertices of the skeleton expansion. 

tioned multiparticle kernels, or, in terms of the standard skeleton expansions, through fully- 
dressed vertices (instead of tree- level ones). For example, in the case of the Q 3 three-gluon 
vertex T afJjU (q, r,p) we have (with a certain hindsight) identified the special leg to be the one 
entering into graph (a x ) of Fig. [T]from the right, carrying momentum q; the corresponding 
vertex SDE is shown in Fig. |2j 

Now, when the Schwinger mechanism is turned off, the various multiparticle kernels 
appearing in the SDE for the Q 3 vertex have a complicated skeleton expansion, but their 
common characteristic is that they are one-particle irreducible with respect to cuts in the 
direction of the momentum q. Thus, for example, diagram (a) of Fig. [3] is explicitly excluded 
from the (gray) four-gluon kernel (b), and the same is true for all other kernels. 

When the Schwinger mechanism is turned on, the structure of the kernels is modified by 
the presence of the composite massless excitations, described by a propagator of the type 
i/q 2 . For example, as shown in Fig. [3j the gray four-gluon kernel is converted into a black 
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FIG. 4: The SDEs for the full three and four gluon vertices in the presence of their pole parts. 
Note that (i) the new SDE kernels are modified with respect to those appearing in Fig. [2j (ii) the 
last term in these SDEs corresponds to the U part of the pole vertices. 

kernel, diagram (&'), which is the sum of two parts: (i) the term (b^), which corresponds to 
a kernel (gray striated) that is "regular" with respect to the g-channel, and (ii) the term 
(b' 2 ), which describes the exchange of the composite massless excitation between two gluons 
in the g-channel. 

Thus, when the replacements of Eq. ( 12 . 1 2 j) and Eq. (12. 8p are carried out, the SDE for the 
different interaction vertices in the presence of their pole parts will be given by expansions 
such as those shown in Fig. HI 

These modifications in the composition of the kernels give rise precisely to the vertices 
V mentioned earlier. A closer look at the structure of the terms comprising the last term in 
Fig. H] reveals that the Lorentz index a (of the leg carrying the momentum q) is saturated 
precisely by the momentum q. Similarly, Bose symmetry forces the same behavior on the 
other two channels, so that, in the end, we obtain a totally longitudinal structure, i.e., 
Eq. (EH]). 

At this point it is advantageous to make the nonperturbative pole manifest, and cast the 
last term of the three-gluon vertex in the form of Fig. HJ by setting 

U aflu (q,r,p) = I a (q)—B^ v (q,r,p) ; (3.1) 
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if abc B(q,r,p) 



ii 




{l,r,p,l) 



FIG. 5: Diagrammatic representation of the transition amplitude I a (q), an d the Feynman rules for 
the different effective vertices -B. 

J a (g) denotes the transition amplitude that mixes a quantum gluon with the massless exci- 
tation, i/q 2 corresponds to the propagator of the massless excitation, and B is an effective 
vertex describing the interaction between the massless excitation and gluons and/or ghosts. 
In the standard language used in bound-state physics, B represents the "bound-state wave 
function" (or "Bethe-Salpeter wave function"). Clearly, due to Lorentz invariance, 



and the scalar cofactor, to be referred as the "transition function" , is simply given by 



Furthermore, notice that the transition amplitude I a (q) is universal, in the sense that it 
constitutes a common ingredient of all V vertices, namely 




(3.2) 




(3.3) 




(3.4) 



14 



• 

(A) vwv^T 
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+ 




(C) 
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FIG. 6: Basic steps in the derivation of the BSEs that govern the dynamics of the effective vertices 
D. 

thus, the difference between the various U a {...} vertices is solely encoded into the structure 
of the wave functions Bs\. The diagrammatic representation of the transition amplitude 
as well as the different effective vertices B are shown in Fig. [5l It is important to recognize 
that the U a {...} depend quadratically on the since I a (q) depends linearly on them. 

As has been explained in the literature, the dynamics of the B functions may be deter- 
mined, at least in principle, from a set of homogeneous (coupled and nonlinear) integral 
equations, known as BSEs. This particular set of equations must admit nontrivial solutions, 
which, when properly adapted to the kinematic details of the problem at hand, will furnish 
the momentum dependence of the wave functions B. 

The way to obtain this set of equations is by first rearranging the vertex SDE in such 
a way as to replace the bare vertex appearing in the original expansion (e.g., Fig. [2]) by a 
fully dressed one (Fig. [6l line A). At the same time, and as a consequence of this rearrange- 
ment, one must replace the standard SD kernel by the corresponding Bethe-Salpeter kernel 
(denoted by /C in Fig. [6]); the two kernels are diagrammatically different, and are formally 



related by a standard all-order formula 58M60| . The next step is to separate the full vertex 



into the "regular" part, namely the part that behaves as a regular function in the limit 
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q — > 0, and the pole part, l/q 2 , as shown in Fig. [6j line B. Note that the full regular part is 
the sum of T m and the term 1Z coming from V, since neither of these two terms diverges in 
the aforementioned limit. Then, this separation is carried out on both sides of the equation; 
since the vertex on the right-hand side (rhs) is now fully dressed, it too possesses a pole l/q 2 . 
The final dynamical equation for B emerges by equating the coefficients multiplying the pole 
term on both sides, as seen in Fig. [6j line C. In the same way, the dynamical equation for 
the regular part is given by the remaining terms. 



IV. RELATING THE GLUON MASS WITH THE TRANSITION AMPLITUDE 

The aim of this section is to derive the fundamental formula that relates the effective 
gluon mass with the square of the transition amplitude. 

To that end, let us go back to the SDE for the gluon propagator with the replacements 
given in Eqs. (12 . 8 j) - (12 . 1 21) already implemented, namely (Landau gauge) 

5 

[q 2 J m (q 2 ) - m 2 {q 2 )]P, v {q) = <?PM + * $>*U ( 4 - X ) 

i=i 

It turns out that the most expeditious way for deriving the gluon mass equation, and from 
it the desired relation between m 2 (q 2 ) and I{q 2 ), is to identify, on both sides of Eq. (14.11) . 
the cof actors of the tensorial structure q^qu/q 2 that survive the limit q 2 — > 0, and then set 
them equal to each other. In doing so, it is clear that the left-hand side (lhs) of Eq. ( 14.ip 
furnishes simply 

[lhsU = q -^m 2 (q 2 ). (4.2) 

On the other hand, the corresponding contribution from the rhs is directly related to the U 
part of the V vertices (see discussion in subsection IIIBj) . which, due to their very definition 
[see Eq. ( I3.4p ], are all proportional to q v . 

The procedure described above is exemplified in Fig. [7] for the particular case of the 
diagram (a[). Specifically, in the first step one separates the regular part T m from the pole 
part V of the full Q 3 vertex T'. In the second step the pole part V is written as the sum of 
the U part (containing the explicit g-channel massless excitation) and the 1Z part. Finally, 
due to the special structure of the IA part, this contribution is proportional to q^q u /q 2 , and 
contributes to the rhs of the mass equation. 
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FIG. 7: Procedure for isolating the contribution of diagram (a^) to the gluon mass equation. 



The next step is to carry out this procedure to the U parts of all diagrams, and determine 
the complete contribution to the rhs of the mass equation, as shown pictorially in Fig. |SJ It 
is clear that from all diagrams containing the U parts (first line of Fig. [H]) one may factor 
out the common quantity I(q 2 ), since, as we have emphasized in the previous section, I(q 2 ) 
is universal (second line of Fig. [8]). Then, quite interestingly, the sum of the terms in the 
parenthesis is nothing else than the diagrammatic representation of I(q 2 ), given in Fig. [5] 
(note that all combinatorial factors work out exactly). 

So, applying this procedure for each one of the fully-dressed vertices appearing in the 
SDE Eq. (14. ip . one can put together the contributions of the several U parts to the gluon 
self-energy, as shown in Fig. |8l 

5 



D 

t=i 



q~ 



(4.3) 



Then, using Eq. f 13 . 2 1) . together with the fact that I u (—q) = —I v (q), one obtains the following 
result 



(4.4) 



Therefore, equating Eq. (14. 2 p with Eq. ( 14. 4p we find that the effective gluon mass is related 
to the transition amplitude through the simple formula (Minkowski space) 



mV) = *rJV)- 



(4.5) 



This last formula may be passed to Euclidean space, using the standard conversion rules, 
together with the expression for I(q 2 ) given in Eq. (I8.17p . In particular, the transition to the 
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FIG. 8: General procedure for determining the contributions of the rhs of Eq. (I4.ip to the effective 
gluon mass. Note that the vertices appearing in this figure correspond to the U parts of the full 
vertices V. 

Euclidean space proceeds by using the standard formulas that allow the conversion of the 
various Green's functions from the Minkowski momentum q 2 to the Euclidean q 2 = —q 2 > 0; 
specifically 



A E (ql) = -A(-(£); m 2 E (q 2 E )=m 2 (-q 2 E ); G E (q 2 E ) = G(-q 2 E ); f = if . (4.6) 

J k J kg 

As a consequence, we have I e (q e ) = ~ H~ Qe)i so ^ na ^ 

m 2 E (q 2 E )=g 2 I 2 E (ql). (4.7) 

An immediate important implication of this last relation is that the gluon mass obtained is 
a positive-definite function for all values of the Euclidean momenta, as expected on physical 
grounds. 

Let us finally mention that the transversality of the full gluon self-energy [see Eq. (12.111) 
guarantees that if one were to consider the part of the mass equation proportional to g^, 
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one would eventually obtain exactly the same relation given in Eq. ( 14. 5p . Note, however, 
that the corresponding derivation is far more subtle and laborious, and hinges crucially on 
the judicious use of a special integral identity (for a detailed treatment of this issue, see [21]). 



V. THE BQI OF THE TRANSITION AMPLITUDES: SDE DERIVATION 

We will now repeat the construction of the last section using the SDEs of the PT-BFM 
formalism. Specifically, we will derive the relations analogous to Eq. (I4.5p . which, in con- 
junction with the BQIs connecting the conventional and PT-BFM gluon propagators, will 
furnish a nontrivial relation between the corresponding transition amplitudes. 



A. General considerations 



In the PT-BFM formalism the natural separation of the gluonic field into a "quantum" 
(Q) and a "background" (B) gives rise to an extended set of Feynman rules, and leads to an 
increase in the type of possible Green's functions that one may consider. In the case of the 
gluonic two-point function, in addition to the conventional QQ gluon propagator, A, two 
additional quantities appear: the QB propagator, A, mixing one quantum gluon with one 
background gluon, and the BB propagator, A, with two background gluon legs. It turns 
out that these three propagators are related by the all-order identities (referred to as BQIs) 

A(g 2 ) = [1 + G{q 2 )]K{q 2 ) = [1 + G(q 2 )] 2 A(q 2 ). (5.1) 

The function G(q 2 ), whose role in enforcing these crucial relations is instrumental, is 
defined as the g^ v form factor of a special two-point function, given by (see Fig. [9]) 

KM = ~W 2 C A [A-(k)D(q-k)H va (-q,q-k,k) 

Jk 

= ^G(g 2 ) + M^L(g 2 ), (5.2) 

where Ca denotes the Casimir eigenvalue of the adjoint representation (N for SU(N)), 
d = 4 — e is the space-time dimension, and we have introduced the integral measure 

. a' 



k (2«yj ik < < 5 - 3 > 

with n the 't Hooft mass. In addition, D ab (q 2 ) = 5 ab D(q 2 ) is the ghost propagator, and H v 



is the gluon-ghost kernel 
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621 ] . The dressed loop expansion of A and H is shown in Fig. [91 
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FIG. 9: Definitions and conventions of the auxiliary functions A and H. 

Notice that the standard ghost-gluon vertex is obtained from H Ufl simply through the 
contraction 

q u H^(q,p,r) = -Tn(r,q,p). (5.4) 

Finally, in the Landau gauge only, the form factors G(q 2 ) and L(q 2 ) are linked to the ghost 
dressing function 

F(q 2 ) = q 2 D(q 2 ) (5.5) 



by means of the all-order relation 571. 163H65] 



F-\q 2 ) = l + G(q 2 ) + L(q 2 ), (5.6) 

which will be extensively used in the ensuing analysis. 

Returning to Eq. (15. ip . it is important to recognize that the two basic functions consti- 



tuting the gluon propagators, namely J m (g 2 )and m 2 (q 2 ), satisfy Eq. (15.11) individually |l3 |. 
In particular, the corresponding masses are related by 

fh 2 (q 2 ) = [1 + G(q 2 )}m 2 (q 2 ) = [1 + G(q 2 )} 2 m 2 (q 2 ). (5.7) 

Finally, in order to obtain from the SDEs of the PT-BFM propagators a mass formula 
analogous to that of Eq. (14. 5p . one needs to introduce the appropriate transition amplitude 
connecting the background gluon with the massless excitation. This new transition ampli- 
tude, whose diagrammatic representation is shown in Fig. [TOl allows us to write the U part 
of the corresponding pole vertices in the BFM as 

U a{ ... } (q,...) = I a (q)^B { ... } (q,...). (5.8) 
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FIG. 10: Diagrammatic representation of the background transition amplitude I a {q)- 

Observe that only the transition amplitude is modified in this expression with respect to 
Eq. ( 13. 4 p when we go to the BFM. This is so because the only background field is the one 
carrying the momentum q, while all other fields are quantum (i.e., they are common to the 
PT-BFM and conventional descriptions). 

B. Relating the transition amplitudes through the SDEs of the PT-BFM 

Consider the SDE for the QB propagator, connecting a quantum and a background field, 
given by the expression 

6 

A-\q 2 )P, v {q) = q 2 P» v {q)+iY<^)^ ( 5 ' 9 ) 

i=i 

with the diagrams (a*) shown in Fig. [TTJ Note that the diagrams are separated into 
three blocks, each of which is individually transverse; this special property of "block-wise" 
transversality is particular to the PT-BFM scheme. 

Let us now apply the diagrammatic procedure shown in Fig. [5] to the new set of diagrams 
appearing in this SDE. Evidently, the quantity to be factored out from all diagrams com- 
prising the W-related part of the SDE is I{q 2 )- Then, it is easy to recognize that the (sub) 
diagrams composing the other factor coincide precisely with those shown in the first line of 
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FIG. 11: SDE obeyed by the gluon self-energy with one quantum (Q) and one background gluon 
(B) entering. Each box encloses a set of diagrams forming a transverse subgroup. The small 
gray circles appearing on the external legs (entering from the right, only!) are used to indicate 
background gluons. 

Fig. [5j thus giving rise to I{q 2 ) again. Thus we arrive at the following expression for the 
mass of the QB propagator, 

m 2 (q 2 )=g 2 I(q 2 )I(q 2 ). (5.10) 

At this point one may use the BQI of Eq. (15. 7\\ to replace m 2 (q 2 ) in favor of m 2 (q 2 ) on 
the lhs of Eq. (15.101) , namely 

Note that this last substitution is legitimate, because the corresponding SDEs have been con- 
sidered in their full, untruncated version (all diagrams kept); therefore the masses obtained 
from them are the same exact quantities that satisfy the BQI. 

Then, direct comparison with Eq. f 14 . 5 1) furnishes the central relation 



I{q 2 ) = [1 + G{q 2 )]I{q 2 



(5.12) 



Interestingly enough, this relation emerges as a self-consistency requirement between the 
results obtained from two formally different, but physically equivalent, versions of the gluon 
propagator SDE. 

Furthermore, due to the special structure of Eq. (15. 8p , where only the transition amplitude 
knows if the g-leg is quantum or background, the result Eq. (15. 12[) implies that the part U 
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FIG. 12: SDE obeyed by the gluon self-energy with two background gluons (B) entering. The 
graphs inside each box form a transverse subgroup. External background legs are indicated by the 
small gray circles. 

of the pole vertices must also satisfy the same BQI, namely (suppressing indices) 

U = [1 + G(q 2 )]U. (5.13) 

Finally, let us consider for completeness the SDE for the BB propagator connecting two 
background gluons, given by the expression 

10 



A- 1 (q 2 )P, u ( q )= q 2 P, u (q)+ l J2( 



"l J fJ,V1 



(5.14) 



i=i 



where the diagrams (a^) appearing on the rhs of the equation are shown in Fig. [121 arranged 
again into individually transverse blocks. 

Repeating the procedure of Fig. E] for this last SDE, we see that, as in the QB case, 
the common quantity to be factored out from all diagrams is again I(q 2 )- Then, the corre- 
sponding sum of (sub) diagrams coincides precisely with those of Fig. [Hll thus giving rise to 
another I(q 2 ). As a result, one obtains 



™V) = <? 2 /V)- 



(5.15) 
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It should be easy to verify at this point that the direct use of Eq. (15. 7p and subsequent 
comparison with Eq. (14. 5 p [or with Eq. (I5.10P ], furnishes again the result of Eq. (I5.12p . 



VI. THREE-GLUON POLE VERTEX 

In the previous section the basic relations Eq. (I5.12p and Eq. (15.131) have been derived 
at the level of the SDEs, by appealing to the general properties of the massless bound state 
formalism and the BQIs relating the gluon propagators in the PT-BFM formalism. In the 
next two subsections we will derive the same relations from the pole parts of the Q 3 vertex 
and the BQ 2 vertex. Specifically, in the first subsection we derive in detail the closed form of 
these vertices, while in the next we identify their ti parts, and carry out a direct comparison. 



A. Explicit construction 



As has been explained in detail in the early 16l420l] and recent literature [2l|, |66] , in order 
to preserve the gauge invariance of the theory in the presence of gauge boson masses, the 
usual vertices must be supplemented by pole parts. In particular, the pole part of the BQ 2 
vertex must satisfy the following WI and STIs 
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66fl, 



q a V afM/ (q,r,p) = m 2 (r 2 )P^(r) - m 2 (p 2 )P^ u (p), 



r»V afiU (q,r,p) = F{r 2 )[m 2 (p 2 )P^p)H^(p,r,q) - m 2 (q 2 )P^(q)H lll/ (q,r,p)}, 
p"V aia ,(q,r,p) = F(p 2 )[m 2 (q 2 )P»(q)H^(q,p,r) - m 2 (r 2 )P»(r)H ua (r,p,q)}. (6.1) 

The quantity H is given by the same diagrammatic representation as that of H, shown in 
Fig. but with the incoming gluon field replaced by a background one. 

In the pole part of the Q 3 vertex, the leg associated with q a is quantum instead of 
background, and, as a consequence, the Abelian-like WI [first in Eq. (16.11) ] is replaced by an 
STI, namely 



q a V a »„(q,r,p) = F(q 2 )lm 2 (r 2 )P«(r)H au (r,q,p)-m 2 (p 2 )P:(p)H a ,(p,q,r)}. 



(6.2) 



The STIs with respect to the other two legs are those of Eq. (16.11) . but with the "tilded" 
quantities replaced by conventional ones. 
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It turns out that the explicit closed form of the two pole vertices in question, namely V 
and V, may be determined from the STIs they satisfy, and the requirement of complete lon- 
gitudinally, i.e., the condition Eq. ( 12. 91) . Specifically, opening up the transverse projectors 
in Eq. ( 12. 9ft . one can write the entire vertex in terms of its own divergences, 

- ^rWp. + %^ (6.3) 

Note that the last term will not contribute because if we apply the STI's, 

?Vp ff Vfer,p)=0. (6.4) 

So, using Eq. (16. ip to evaluate the various terms, and after a straightforward rearrangement, 
we obtain the following expression for the pole part of the BQ 2 vertex, 

V a , v (q,r,p) = ^[m 2 (r 2 )-m 2 (p 2 )}P^r)P pu (p) 

+ D(r 2 )[m 2 (p 2 )PMH pa (p,r,q) -m^P^P^i^H^q^p)]^ 

+ D(p 2 )[m 2 (q 2 )P^q)H pil (q,p,r) - m 2 {r 2 )Pj:{r)H pa {r^q)}p v . (6.5) 

Applying the same procedure, but using now the STIs of Eq. (I6.2p . together with the 
condition of Eq. (12. 9p . we derive the closed expression for the pole part of the Q 3 vertex, 

Vap,(q,r,p) = D(q 2 )[m 2 (r 2 )H pa (r,q,p) - m 2 (p 2 )H ap (p,q,r)]PP(r)P°(p)q a 

+ D(r 2 )[m 2 (p 2 )P^p)H pa (p,r jq ) - m 2 (q 2 )P^q)P^p)H pa (q,r,p)}r fl 

+ D(p 2 )[m 2 (q 2 )PZ(q)H pix (q,p,r) - m 2 (r 2 )P£(r)H pa (r,p,q)}p„. (6.6) 

Let us discuss next certain issues related to the self-consistency of the previous vertex 
construction. Notice that, in order to obtain the expressions in Eq. ( 16. 5p and Eq. ( 16.61) . 
one must apply sequentially the WI and the STIs. In doing so, the Bose symmetry of both 
vertices is no longer explicit, and the result obtained is not manifestly symmetric under 
the exchange of the quantum gluon legs. Furthermore, seemingly different expressions are 
obtained, depending on which of the two momenta acts first on V or V. However, if one 
imposes the simple requirement of algebraic commutativity between the WI and the STIs 
satisfied by the three-gluon vertex, the Bose symmetry becomes manifest. For example, 
using Eq. ( 16. 5p we can see that the elementary requirement 

gW QM ,(g,r,p) = rYV(?.^), (6-7) 
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gives raise to the following identity 

F(r 2 )P»{p)q a H^{p,r,q) = -r.P^p). (6.8) 

A similar identity is obtained by imposing the requirement of Eq. (16. 7p at the level of V, 
namely 

F(r 2 )Pi:(p)q a H m (p,r,q) = -F{q 2 )P^{p)r a H m {p, q, r). (6.9) 

Quite interestingly, the identities in Eq. (16. 8ft and Eq. (16. 9 h are a direct consequence of 
the WI and the STI that the kernels H and H satisfy, when they are contracted with the 
momentum of the background or quantum gluon leg, namely j^] 

q a H, a (p,r,q) = -^F~V) " V^V), 

q a H, a (p,r,q) = -F(q 2 )\p p F-\p 2 )C(p, qi r) + r a F-\r 2 )H pa (p,q,r)}, (6.10) 

where C(p, q, r) is an auxiliary ghost function. 

Indeed, use of Eq. (16.101) into Eq. (16. 8ft and Eq. (16. 9ft . respectively, leads to a trivial 
identity. Conversely, one may actually derive Eq. (16.101) from Eq. ( 16 .8p and Eq. (\Q.9\i ; for 



example, starting with Eq. (16. 8p . and using also the identities [67 1 



•pPH^ip^q) = r a F 1 (r 2 ) -T a (q,p,r), 

q a f Q (q J p J r)=p 2 F- 1 (p 2 )-r 2 F~ 1 (r 2 ) J (6.11) 

one can easily reproduce Eq. (16.101) . 

Evidently, these constraints allow us to cast the pole part of the BQ 2 vertex into a 
manifestly Bose symmetric form with respect to the quantum legs, 

V a ^{q,r,p) = -|[m 2 (r 2 ) - m 2 (p 2 )]PP(r)P pu (p) + S afiv (q,r,p) - S aup (q,p,r), (6.12) 

with 



S ai u,(q,r,p) = D(r )m (p )PP(p)H pa (p,r, q)r 



_ ^D(r 2 )m 2 (q 2 )P^qM + P^p)]H pa (q,r,p). (6.13) 
Finally, for the pole part of the Q 3 vertex, the Bose symmetric expression reads 

V a ^(q,r,p) = S apu (q,r,p) -S pau (r,q,p) - S Vfla (p,r,q), (6.14) 
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with 



Sanu(q,r,p) 



q y(q 2 )m 2 (r 2 )PP(rM + P°(p)]H pr7 (r, q ,p) 
- f i%VV)^W + P°{r)]H pa {p, q, r). 



(6.15) 



Note the fact that, as it should be, setting in Eqs. (16.12n - fl6.14j) the tree-level values F = 1 
and H pv = g pv for the ghost dressing function and the gluon-ghost kernel, respectively, one 
recovers the expression for the "abelianized" three-gluon vertex first presented in 68] . 



B. Transition BQI derived from the pole vertices 

We next identify from the explicit expressions for the BQ 2 and Q 3 pole vertices the 
terms U and hi. To that end, and in complete accordance with the discussion presented 
in subsection III Bl we apply a kinematic and a dynamical criterion, namely (i) determine 
U by collecting all terms containing the tensorial structure q a /q 2 , and (ii) extract U by 
discarding (and reassigning them to W) all terms that, when inserted into the SDE of the 
gluon self-energy, do not survive the q — > limit, i.e., do not contribute to the gluon mass 
equation. The final objective is to infer the validity of the BQI of Eq. (15.131) . connecting U 
and U. 

Applying criterion (i), and denoting by U and U the resulting expressions, it is straight- 
forward to obtain from Eq. (16.121) 

U ap „(q,r,p) = q ^m 2 {r 2 )P^r)P pu {p) - r Jt D {r 2 )m 2 {q 2 M + P>(p)]T p (p, q, r)\ 
r <H- p\ 

' 5 (6.16) 



where Eq. (15. 4p has been applied in order to relate the contractions of the auxiliary ghost 
functions PI with the conventional gluon-ghost vertices T p , and ( r ^J is obtained from the 
term that is shown explicitly after carrying out the exchanges indicated. In addition, for 
the R part of the vertex we find 

Ra^{q,r,p) = r ^F(r 2 ){m 2 (p 2 )PZ(p)H pa (p,r,q) - \m 2 {q 2 )[g p v + P v p (p)]H ap (q,r,p) 



r <H> p 



(6.17) 
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Similarly, from Eq. (I6.14p . and using the special STI Eq. (I6.10p . we obtain 
U a ^(q,r,p) = ^^F(q 2 )m 2 (r 2 )H p(7 (r,q,p)P^r)P^p) 

- ^(r 2 )m 2 (g 2 )^ + P^p)]r p (p,g,r)} - (j^), (6.18) 



and 



R a ^(q,r,p) = r ^F(r 2 )^m 2 (p 2 )PP(p)H pa (p,r,q) - \m 2 {q 2 )\gl + P^p)}H ap (q,r,p) 



r <H> p 



(6.19) 



Let us now assume for a moment that U = U and U = U, as well as R = 1Z and R = 1Z. 
Let us further use Eq. (15 .7R . and substitute m 2 (q 2 ) = m 2 (g 2 )[l + G(q 2 )]~ l into Eq. (I6.18p . 
Then, employing Eq. (I5.6P valid in the Landau gauge, dropping L(q 2 ) since L(0) = [thus 
applying effectively criterion (ii)], it is relatively easy to establish that the (would-be) U 
and U fail to satisfy Eq. (I5.13P due to the presence of a very particular term. Specifically, if 
H pa — > g pa in the first line of Eq. ( 16. 18}) . then the BQI of Eq. ( 15.13P would be satisfied. 

Therefore, write the H pa (r,q,p) in Eq. (16. 181) as 

H P a(r,q,p) = g pa + HQ(r,q,p), (6.20) 

where the g pa corresponds to the tree- level value of H pa , and H® a contains the all-order 
quantum corrections; an exactly analogous expression holds for H ap (p,q,r), with p ■<->• r. 
Note that in both cases the momentum q is carried by the ghost leg, see Fig. [13j Then, if 
we define the part U as [see Eq. (I2.19P ] 

U a ^(q,r,p) = U aiiu (q } r,p) -U^ u (q,r,p), (6.21) 

with 

KM,r,p) = %F{q 2 )m 2 {r 2 )P?{r)P:{p)H%{r,q,p) - (* ** P ) , (6.22) 
it is clear that 

U a ^{q,r,p) = [l + G(q 2 )]U ailu {q,r,p). (6.23) 

The real justification for discarding U^ u from U apbV is provided precisely by invoking 
criterion [ii): the discarded term does not survive the limit q — > when inserted into the 
corresponding gluon SDE. 
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V 

FIG. 13: Diagrammatic representation of Eq. (16.240 . 

To demonstrate that this is indeed so, note that in the Landau gauge 
(k + q) x A\ a (k) = q x A\ a (k), and the momentum of the incoming ghost leg factorizes out 
of the loop integrals. This observation allows us to cast H® in the form 

H%(r,q,p) = q x T Xpa (q,r,p). (6.24) 

where T simply represents the rest of the diagram (see Fig. [T3|) . If we assume now that T 
has a finite (non-divergent) value in the limit q — > 0, then 

KmH%r( r ,<l,p) = limg A T Apff (g,r,p) = 0. (6.25) 
So, using Eq. (I6.24p . the term given in Eq. (16 .22 j) becomes 

ZCfar.p) = ^F(q 2 )m 2 (r 2 )PP(r)P°(p)T Xp(7 (q,r,p) - ( T ^ P ). (6.26) 

When U' aiiV is written in this form, it is clear that its pole 1/q 2 is actually compensated by 
q a Q X , while the remaining terms cancel themselves when r = —p. Therefore, the entire term 
U'apv vanishes as q —> 0, and should not be included in the U part of the Q 3 pole vertex. 

Finally, returning to the identification U = U, note that the expression in Eq. ( 16. 16[) 
does not get modified by the application of criterion (m), because the above argument of 
discarding does not apply in this case. The reason for that is simply the channeling 
of the momenta in H pa (q,r,p) and H pa (q,p,r) (encoded into the gluon-ghost vertices T p ) 
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is different from that of H pcr (r,q,p); the momentum entering in the ghost leg is no longer 
q, but rather r or p, and the corresponding H®s do not satisfy Eq. (I6.25p . As a result, no 
terms need be discarded from U, and, therefore, W = 0. 

Let us now cast U afiu (q,r,p) and U afiu (q,r,p) into their canonical form of Eq. (I3.4p and 
Eq. (15.81) . respectively, namely 

U aiiV {q,r,p) = -^I^B^q.r.p), 

U af ,„{q,r,p) = -^T^B^fa^p); (6.27) 

the minus sign comes from the extra imaginary factor appearing in the Feynman rule of 
the effective vertex, see Fig. |5j The important point to recognize is that B piJ (q,r,p) is 
common to both vertices, because the legs carrying r and p are quantum ones, and the 
difference induced due to the quantum or background nature of the leg carrying momentum 
q is entirely encoded into the form of the corresponding transition amplitudes, I{q 2 ) and 
I(q 2 ), respectively. Then, Eq. ( 16 . 2 7|) and Eq. (I6.23|) imply directly the validity of Eq. (I5.12p . 



VII. DIAGRAMMATIC DEMONSTRATION OF THE BQI 

It is well-known that the BQIs relating the conventional and BFM Green's functions 
are formally obtained by resorting to the powerful BV formalism. On the other hand, 



the PT (or its generalized version 3JJ) furnishes an equivalent diagrammatic derivation, 
which makes extensive use of the STIs satisfied by the kernels appearing in the SDEs of the 
Green's functions in question. In the previous sections the BQI of Eq. (" 15 . 1 2[) relating the 
transition amplitudes I(q 2 ) and I(q 2 ) has been obtained as a self-consistency requirement 
between two SDEs, and from the corresponding BQI relating U and 1A (sections |V] and 
|VT| respectively). The objective of this section is to carry out a PT-guided demonstration 
of this particular BQI, through the systematic conversion of the set of Feynman diagrams 
defining one transition amplitude into the set defining the other (shown in Figs. [5] and 
[TU]) . Essentially this conversion proceeds by allowing the longitudinal (pinching) momenta 
contained in the tree-level three-gluon vertex [graph (di)] to act on the adjacent kernel; 
therefore, as we will see, a crucial ingredient for completing this construction is the knowledge 
of the STI satisfied by the effective vertex B pv . In addition, the requirement that the BQI 
be diagrammatically exact imposes a strong constraint on the set of ghost diagrams that 
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FIG. 14: Detailed representation of diagrams (d±) and (dj). Diagram (d\) has a symmetry factor 
1/2; diagrams (ds) and (^7) vanish in the Landau gauge. 

contribute to the transition amplitude, which may be translated (under mild assumptions) 
into the vanishing of the corresponding subset of B vertices. 

A. The PT construction at the level of the transition amplitude 

Let us start by considering the diagram (di) of the background transition amplitude, 
shown in Fig. [HJ whose contribution is given by 

(di)a = \Ca Jr^ u A^(k + q )A^(k)B pa . (7.1) 

In order to avoid notational clutter we will suppress the arguments of the momenta in the 
vertices, which can be easily recovered from the figures. In what follows, the coupling and 
color dependence of the diagrams comprising the evaluation of I a (q) is restored by setting 
= gS^c 

We know that the tree-level vertex rSjL contains terms which are proportional to so 
that one cannot take directly the limit £ = to work in the Landau gauge. Specifically, the 
tree-level part of the BQ 2 vertex is given by 

FSUg, -* -<?,*) = rgUg, —k -q,k)- ir^fa -* -g,k), (7.2) 

where the purely longitudinal "pinch part" T p reads 

r a^(?' ~ k -1>k) = g ai iK + g au (k + q)^ (7.3) 
31 



and 

rgl(g, -k -q,k) = -(2k + q) a g pv + (k - q) p g av + (2q + k) v g ap (7.4) 

is the tree-level value of the conventional three-gluon vertex. Using then the decomposition 
Eq. (17. 2p . we obtain from Eq. (17. ip 

(dt) a = (dj) - ^C A J T^A^k + qW P (k)B pa . (7.5) 

Now, after the shifts k + q i— » k and k \- > —k, and applying the antisymmetry property of 
the effective vertex, namely B pa = —B ap , the two contributions coming from the pinch part 
of the vertex in Eq. (17.51) sum up and cancel the symmetry factor, giving the result 

(di) a = (di) a - % -C A J k v A up {k)A° a {k + q)B pa . (7.6) 

Therefore, employing the identity 

1 k p 

f^ Vp {k) = ¥ , (7.7) 

the £ _1 term cancels, and we can project Eq. (I7.6P to the Landau gauge (£ = 0), 

~ f k p 

[di) a = (di) a ~ iC A J ~^K{k + q)B pa , (7.8) 

where the gluon propagator assumes a totally transverse form, i.e., A MJ ,(/c) = A(k 2 )P pu (k). 

At this point, in order to evaluate the integral on the rhs of Eq. (I7.8I) . the knowledge of 
the STI satisfied by the effective vertex B pa is required. Observe then that, if we compare 
Eq. (16.16P with Eq. ( 16.271) . the following structure can be identified as the B pu , 



I(q 2 )B pu (q,r,p) = -m 2 \r 2 )P?(r)P pv (p) + !j-D(r 2 )m 2 (q 2 M + P u p (p)]T p (p, q, r) 

r <H> p" 



(7.9) 



Therefore, contracting Eq. (17.91) with respect to the momentum r M , and breaking the trans- 
verse projector P p {p), we find 

I(q 2 )r p B pu (q,r,p) = F(r 2 )m 2 (q 2 )T u (p,q,r) 

- ^m 2 (q 2 )[F(r 2 )pPT p (p,q,r) + F(p 2 yT p (r,q,p)]. (7.10) 

Now, the second term on the rhs does not contribute, since p v — » ik + q) a is annihilated by 
the transverse projector of A^(/c + q) in Eq. (I7.8p . 

32 



Thus, employing the STI Eq. (I7.10p . adapted to our kinematical configuration, i.e., 
r = —k and p = k + q, together with Eq. (I5.10p . Eq. (17. 8p becomes 

{di) a = (d 1 ) a + tg 2 C A I(q 2 ) [ A°(k + q)D(k)T a 

Jk 

= {dx) a - ig 2 C A I(q 2 )q X ! A° a (k + q)D(k)H Xa , (7.11) 

Jk 

where, in the second line, Eq. (15 .4p has been used. 

At this point, one recognizes on the rhs of Eq. (17. lip the appearance of the auxiliary 
two-point function A, defined in Eq. (15.21) ; so, we obtain 

(di) a = {d l ) a + I{q 2 )q x Kx{q) 

= (d 1 ) a + [G(q 2 ) + L(q 2 )]I a (q). (7.12) 

Finally, since L(0) = 0, one can drop the term L(q 2 ), exactly as was done in the previous 
section (and with the same justification), finally arriving at 

{d 1 ) a = {d l ) a + G{q 2 )I a {q). (7.13) 

Let us next observe that (i) diagrams (ofa) and (di) in Fig. [5] can be converted automati- 
cally to background ones, given that the tree-level vertex BQ 3 is identical to the conventional 
four gluon vertex Q A , and that (ii) diagrams (d 6 ) and (0^7) vanish in the Landau gauge, see 
Fig. dU this happens because, after integration, the contribution of the ghost loop nested 
inside them will be proportional to k p , and it will vanish when contracted with the A pp /(fc). 
So, 

(d 3 ) = (d 3 ), A) = ((k), (4) = (d 7 ) = 0. (7.14) 

Therefore, comparing Eq. (I7.13P and Eq. (17.141) with Eq. (I5.12p . we see that, in order for 
the full BQI to be diagrammatically realized, the relation 

(d 2 ) = (d 2 ) + (d 5 ) + {ds) (7.15) 

should be satisfied. 

B. A constraint on the ghost sector 

In the previous subsection we have derived Eq. (17.151) from the requirement that the basic 
BQI Eq. (I5.12p be diagrammatically realized. The common characteristic of all diagrams 
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l + k 




FIG. 15: Detailed representation of the diagrams (^2)5 (^5), an d ids), appearing in Eq. (|7.15p . 

comprising Eq. fl 7 . 1 5 [) is that they contain vertices B which mix the nonperturbative massless 
bound-state with two ghosts. Given that the nature of these vertices is practically unknown, 
it would be interesting to inquire what restriction Eq. ( I7.15P may impose on their structure, 
and whether this restriction is compatible with other pieces of nonperturbative information 
on the ghost sector of Yang-Mills theories. 

In the formulas that follow, color factors are shown explicitly, while the arguments of 
the momenta in the vertices are suppressed as before. Then, the diagrams appearing in 
Eq. (17.15P are given by 



{d 2 ) ab = ~igC A 5 ab I T^D(k)D(k + q)B 
Jk 

{d 2 )t = ~igC A 6 ab [ f^D(k)D(k + q)B 
Jk 



(d b ) ab = g 3 C 2 A 5 ab [ T^D(k)D(k + q)B ! D{k + l)A^(l)T v 

Jk Jl 



(dsT a = g 2 f adx r m / D(k + q)D(k + l)^ a il)B b u cmd . (7.16) 

Jk J I 

Observe that the pole part of the internal gluon-ghost vertex T v does not contribute to (efe) 
in the Landau gauge, due to the longitudinality condition that it ought to satisfy. Therefore, 
only the regular part of this vertex survives in (ds). 



We will next use the SDE satisfied by the ghost propagator D(k), written as 35] 

g 2 C A D(k) jD{k + l)^{l)Y v = ik a [D(k) - D^(k)), (7.17) 

to cast (fife) in the form 

{d 5 )f = igC A 6 ab [ k a [D{k) - Di°\k)]D{k + q)B. (7.18) 

Jk 
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Substituting the above results into Eq. (I7.15p . we arrive at the integral constraint 

iC A b ah [ ^D(k + q)B = gf°*>f*<*n [ fD(k + q)D(k + I) A u a (l) B h u cmd , (7.19) 

Jk & Jk J I 

relating B and B v . 

We next turn to the full BFM gluon-ghost vertex, T, which satisfies the following all-order 



wi ea, 



q a T a (q,r,p) = D-\r 2 ) - D'^p 2 ). (7.20) 



Let us suppose that, due to nonperturbative dynamics, T acquires a pole part; then, in 
full analogy to Eq. ( 12. 8p . we define 

f' a (g,r,p) = f a (q,r,p) + V a (q,r } p), (7.21) 

as the sum of a regular part and a pole part. In addition, and again in analogy with 
Eq. f)2.12p . we will assume that also the ghost propagator acquires a mass term, to be 
denoted by ml(q 2 ), and so 

D-\q 2 ) = q 2 F-\q 2 ) D7 m \q 2 ) = q 2 F-\q 2 ) - m 2 c (q 2 ). (7.22) 

Then, it is clear that we must have 

q Q T a (q,r,p) = r 2 F m \r 2 ) - p 2 F-\p 2 ) , 

q a V a (q,r,p) = m 2 c (p 2 ) - m 2 c (r 2 ) , (7.23) 

so that finally 

q a T' a (q,r,p) = D m \r 2 ) - D m \p 2 ) , (7.24) 

as it should. 

If we now assume the existence of a massless pole in the g-channel which, in addition, 
satisfies the longitudinality condition 

pP(q)Mq,r,p) = 0, (7.25) 

then, the WI obeyed by V a (q,r, p), may be solved, yielding 

V a (q,r,p) = ^[m 2 c (p 2 )-m 2 c (r 2 )}. (7.26) 

Note that V a (q,r, p) is purely of the U type (no 7Z terms). 
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Then, using the corresponding expression Eq. (15. 8p for this vertex 

V a (q,r,p) = -^I(q 2 )B(q,r,p), (7.27) 
q 2 

and equating with Eq. (17.261) . we obtain a relation between the effective vertex B and the 
mass-like term in the ghost propagator 

I(q 2 )B(q,r,p) = m\{r 2 ) - m 2 c (p 2 ). (7.28) 



Now, large- volume lattice results [4H7J, |9j, together with a variety of analytic studies 
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43|, 



451 ]. establish that (i) the ghost propagator in the Landau gauge remains massless, which 
means that m 2 (q 2 ) = 0. (ii) The dressing function F m (q 2 ) is infrared- finite; this happens 
essentially due to the fact that the gluon mass m 2 (q 2 ) regulates it, in the same way as 
J m , namely through the qualitative replacement lng 2 — > ln(g 2 + m 2 ) [see discussion after 
Eq. f[2T2]) ]. Therefore, since m 2 c {q 2 ) = but I(q 2 ) ^ 0, Eq. flT^HD implies that B = 0. 

Returning to the constraint of Eq. (17. 19[) and setting B = 0, we have that the integral on 
the rhs involving B u must vanish also. Even though from this statement one cannot math- 
ematically conclude that B v must vanish identically, this seems to be the most reasonable 
scenario. Indeed, the rather fine-tuned alternative would invoke a subtle conspiracy between 
the various form factors comprising B u , which ought to combine in such a way (and change 
signs correspondingly) as to make this integral vanish. 



VIII. COMPARISON OF THE TWO MASS -GENERATING FORMALISMS 

In this section we first derive the full expression for the transition amplitude in the Landau 
gauge, and then we show that, when judiciously combined with Eq. (14.51) and Eq. (I8.19p . it 
gives rise to an integral equation for the gluon mass, which coincides exactly with the one 
obtained within the SDE formalism of [2l|] . We then proceed to a brief study of how the 
distinct approximations employed within the two formalism lead to differences in the form 
of the gluon masses obtained. 



A. Exact formal equivalence 

From the analysis presented in the previous section we have concluded that diagrams 
with ghost loops do not contribute to the transition amplitude. Furthermore, one may 
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FIG. 16: Configuration for the diagrams of I a (q) that contain only gluon loops. The symmetry 
factors for the diagrams are S(di, d$, d^j = (1/2, 1/6, 1/2). 

demonstrate that diagram (ofa) in Fig. [161 which contains the effective vertex B pT p mixing 
three quantum gluons with the bound state, vanishes in the Landau gauge. Effectively, with 
the configuration shown in Fig. [161 we can write down the contribution of this diagram, 

= \J J? { $Sr*V{h + q)Kl(i + k)£&q)E%$. (8.1) 

As in the preceding sections, the arguments of the momenta in the vertices have been 
suppressed. Observe that, in order to evaluate this diagram in the Landau gauge, one 
should know how three transverse projectors act over the effective vertex B pT p. Consider 
then the decomposition Eq. (I2.2ip for the BQ 3 pole vertex, 

Vxpr/s(q,r,p,l) = U\prp(q,r,p,l) + 7^ Apr(3 (g, r,p, Z), (8.2) 

where, according with the discussion carry out in Section [TTl the 1Z part of the vertex satisfies 
the transversality condition 

P°P(r)P» T (p)P^(l)n XpT p(q,r,p,l) = 0. (8.3) 

Thus, when Eq. f)2.9p is applied for the case of the BQ 3 pole vertex, namely, four transverse 
projectors canceling the vertex, we obtain from Eq. (18. 2p the following result 

P^(r)P^(p)P^(l)U apT p(q } r,p,l) = q ^P^(r)P^(p)P^(l)q x V XpT0 (q,r }P} l). (8.4) 

On the other hand Eq. ( 15. 8ft becomes for the U part of this pole vertex, 

U ap ri3(q,r,p,l) = ig-^I(q 2 )B pr p(q,r,p,l). (8.5) 
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Therefore, applying three transverse projectors on Eq. (18. 5p and equating the result with 
Eq. ( 18. 4 p we can relate the effective vertex B pT p with the WI satisfied by the BQ 3 pole vertex 
when it is contracted with respect to the momentum of the background gluon leg, i.e., 

igI(q 2 )P^(r)P^(p)P^(l)B pT p(q,r, P J) = P^(r)P^(p)P^ (8.6) 

Once the above connection has been established through Eq. (18. 6p . one may repeat the steps 
presented in subsection VI B of [21], thus demonstrating the vanishing of diagram (d 3 ). 

Therefore, the full transition amplitude will be given in the Landau gauge solely by the 
sum of diagrams (di) and (d^). Thus, applying the conventions of Fig. [161 we obtain for 
diagram (di) the expression, 

(di)a = l -C A jv^ v ^{k + q)A»»{k)B pa . (8.7) 

One observes that the above integral has one free Lorentz index, which only can be saturated 
by the external momentum q. So, using the elementary WI for the tree-level three-gluon 
vertex, 

q X Tfl{q, -k -q,k) = k 2 P, u (k) - (k + qfP^k + q), (8.8) 

and after the appropriate shifts in the integrated momentum, we deduce in the Landau 
gauge the result, 



= %<7 A K)a = %C A % [ k 2 A;(k)A^(k + q)B pa . (8.9) 

q q Jk 



Now, diagram (d 4 ) contains the tree- level four-gluon vertex r^ ^. After the color algebra 
and using the standard Feynman rule for this vertex, we obtain for the prefactor of the 
diagram, 



\gf bcm r n r { Z™ = -/Cl{g av g„ - g a ^)5 ab . (8.10) 



Thus, we get in the Landau gauge the following expression, 

(d 4 ) a = \g 2 C\{g av g^ - g a pg^)J^(k + q)A^(k)Y^(k)B p)J , (8.11) 
where we have defined the loop integral 

Yf(k) = jA"\l)A^{k + l)T yrX . (8.12) 
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Note that in the Landau gauge, the fully-dressed three-gluon vertex appearing in this integral 
only contains the regular part, since the transverse projectors of the gluon propagators 
trigger Eq. (12.91) for that vertex. 

As before, the integral in Eq. (18. lip only can be saturated by the external momentum q. 
Moreover, due to the Bose symmetry of the three-gluon vertex is straightforward to show 
that the integral Eq. (I8.12p is antisymmetric under the v •<-)■ /3 exchange, and given also the 
antisymmetry of the prefactor under the same exchange, one can write 

Y^(k) = (k» g e-keg: / )Y(k 2 ) ; Y = ~^g^ {k). (8.13) 

With these observations, Eq. (18.111) can be cast in the form 

(d 4 )a = V(4)a = ld 2 C 2 A % [ [{kq)g m + q,q,]Y{k 2 )A^{k + q)A^(k)B pa . (8.14) 

Employing now the results obtained for diagrams (di) and (0^4), we derive the complete 
expression of the transition amplitude in the Landau gauge [see Eq. (13. 3p ]. 

I(q 2 ) = ^[( dl ) + (d 4 )] a = ^C A [ k 2 A^k)A^(k + q)B pa 
1 1 Jk 

+ I IWffw + Q^]y(k 2 )^(k + q)A^(k)B pu . (8.15) 

1 Q Jk 

To proceed further, we decompose the effective vertex B pa in the tensor basis 

B pa = B x g pa + B 2 q p q a + B 3 (k + q) p {k + q) a + B 4 k p q a + B 5 k p (k + q) a . (8.16) 
When this decomposition is inserted in Eq. (I8.15p . only the form factor B\ survives, since 
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26j was shown that B 2 = and the rest of form factors are canceled in the Landau gauge 



by the transverse projectors. Therefore, Eq. (I8.15P becomes 

I(q 2 ) = - 2 C A [ k 2 A^k)A^k + q)B 1 
Q Jk 

+ |^ / [{kq)g m + q pqi ]Y{k 2 )A^k + q)A pr/ (k)B l , (8.17) 
A Q Jk 

which, quite remarkably, allows to express the full transition amplitude, for general value 
of q 2 , in terms of one single form factor, namely, B\. At this point we arrive at the crucial 
observation which enables us to relate the mass equation obtained in the context of the 
PT-BFM formalism with the bound-state formalism showing that, indeed both formalisms 



39 



m 2 (q 2 ) 



q- 



x 




+ 



V 




FIG. 17: Concise diagrammatic representation of the combined operations leading to the all-order 
gluon mass equation; m 2 (q 2 ) is defined in Eq. ([57 



are self-consistent and interconnected. From Eq. (I7.9|) one may obtain, after identify the g pu 
component, a closed expression for the form factor B\ in terms of the gluon mass 

I(<Z 2 )fli(<7,r,j)) = [l + G^I^B^^p) 

= m 2 {p 2 ) -m 2 (r 2 ), (8.18) 



where the BQI Eq. (15.121) has been applied. Therefore, using this relation and after a 
straightforward rearrangement Eq. (I8.17P yields 

iC A 1 



/ 2 (^ 



1 + G{q 2 ) q 2 



(8.19) 



In this expression, the quantity 



JCZ(q,k) = g^[(k + q) 2 -k 2 }h + -ig 2 C A [Y(k + q)+Y(k)} 
+ -ig 2 C A {q 2 g^ - 2q^q^)[Y(k + q) - Y{k)\. 



(8.20) 



represents the SD kernel obtained from the sum of the two graphs on the rhs of the equation 
depicted in Fig. [T71 Thus, using the mass formula Eq. ( 14. 5 h on the lhs of Eq. ( 18.19!) . we 
arrive to the result 

1 



2/ 2^ 
m [q 



ig 2 C A 



1 + G(q 2 ) q 2 



m 2 (k 2 )A, p (k)AP(k + q)}CZ(q,k) 



(8.21) 



Quite remarkably, this result coincides exactly with the full mass equation derived in 21 ] 
from the SDE of the gluon propagator, following the procedure shown pictorially in Fig. [TTl 
whose formal aspects are considerably different compared to those of the massless bound- 
state formalism. The emergence of an identical dynamical equation for the gluon mass 
provides an impressive self-consistency check between these two formalisms. 
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B. Differences in the numerical implementation 



26) 



At this point it is important to recognize that, although formally equivalent, the two ap- 
proaches ["SDE" vs "massless bound-state"] entail vastly different procedures for obtaining 
the desired quantity, namely the functional form of m 2 (q 2 ). Given that in practice approx- 
imations must be carried out to the fundamental equations of both formalisms, the results 
obtained for m 2 (q 2 ) will be in general different; indeed, the formal coincidence proved above 
is only valid when all equations involved are treated exactly. 

To appreciate this issue in some quantitative detail, recall that, within the massless 
bound-state formalism, the starting point for obtaining the momentum- dependence of the 
gluon mass is Eq. ( I8.18p . Specifically, in the limit q — >■ one obtains the following exact 
relation for the derivative of the effective gluon mass (Euclidean space) 

= -T(0)B[(p 2 ), (8.22) 

where the prime denotes differentiation with respect to (p + q) 2 and subsequently taking the 
limit q — > 0, i.e., 

It turns out that the function B[(p 2 ) satisfies its own homogeneous BSE [see Section HTT1 
Fig. [6] and discussion below] , of the general form 

B[(p 2 ) = [ )C(p,k)B[(k 2 ), (8.24) 
Jk 

where K. corresponds to the Bethe-Salpeter four-gluon kernel, shown in line A of Fig. [TS] 
[see also Fig. [6] . Evidently, the exact treatment of this equation would require the complete 
knowledge of the four-gluon kernel, which, however, is a largely unexplored quantity. There- 
fore, in order to obtain an approximate solution to Eq. (I8.24p . one resorts to the "improved" 
ladder approximation of this kernel, shown diagrammatically in Fig. [THJ by dressing the 
gluon propagators but keeping the vertices at tree-level. Then, the numerical solution of 
Eq. (I8.24p yields for B[(q 2 ) the function shown on the left panel of Fig. [JUJ 

With this information at hand, the gluon mass may be obtained through direct integration 
of Eq. (18T22D . namely 

~ r 2 

m 2 {q 2 ) = m 2 (0) - 1(0) / dxB[(x). (8.25) 

Jo 
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FIG. 18: (A) The kernel /C of the BSE satisfied by B[ and the approximation introduced in 26]. 



(B) The full kernel /C SD of the mass equation, and the approximation employed in 



As for the constant 1(0), using the BQI Eq. (I5.12p . as well as Eq. (I5.6P and Eq. (14. 5p . one 
obtains 



7(0) = F- 1 (0)/(0) =F- 1 (0), 



'm 2 (0) 



'A-^O) 



(8.26) 



g* V 47ra s 

which allow us to estimate 7(0) from the lattice values of the ghost dressing function and 
the gluon propagator at zero momentum, treating a s as an adjustable parameter. Thus, 
one finally arrives at the dynamical gluon mass shown on the right panel of Fig. [19] [red 
(continuous) curve] 



Let us now turn to the SDE approach 



21], and compare with the corresponding ap- 



proximations and numerical results. In this case the basic dynamical equation is that of 
Eq. (I8.2ip . whose central ingredient is the quantity Y(k 2 ), entering into the kernel JC SD , 
given in Eq. (18. 20ft . The quantity Y(k 2 ), shown in line B of Fig. [THJ involves the fully- 
dressed three-gluon vertex, which too is rather poorly known. As a result, in {21] the quan- 
tity Y(k 2 ) has been approximated by its perturbative (one-loop) expression, by replacing 
the fully dressed internal gluon propagators and three-gluon vertex by their tree-level values, 
as shown in Fig. [TBI Under these approximations, the numerical treatment of Eq. (18.21)) 
gives rise to the solution shown on the right panel of Fig. [19] [blue (dotted) curve] . 

Note that both masses have been normalized in such a way as to coincide at the origin; 
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FIG. 19: (Left panel): The general form of B' 1 (q 2 ) obtained from the numerical solution of 
Eq. (|8.24|) . in the improved ladder approximation. (Right panel): The two solutions obtained for 
the gluon mass within the massless bound-state [red (solid)] and SDE [blue (dotted)] approaches 

in particular, m 2 (0) = 0.14 GeV 2 , which is the saturation point reached by the gluon 
propagator simulated on the lattice, when renormalized at /j, = 4.3 GeV (the last point 
available in the ultraviolet tail). It is clear from this direct comparison that the two formally 
equivalent approaches lead to qualitatively similar results, which, however, do not coincide, 
due to the inequivalence of the approximations employed. 



IX. CONCLUSIONS 



In this article we have studied in detail the mechanism of gluon mass generation within the 
massless bound-state formalism, which constitutes the formal framework for the systematic 
implementation of the Schwinger mechanism at the level of four- dimensional non-Abelian 
gauge theories. The main ingredient of this formalism is the dynamical formation of massless 
bound-states, which give rise to effective vertices containing massless poles; these latter 
vertices trigger the Schwinger mechanism, and allow for the gauge-invariant generation of 
an effective gluon mass. The principal advantage of this approach is its ability to relate the 
gluon mass directly to quantities that are usually employed in the physics of bound-states, 
such as transition amplitudes and bound-state wave functions, as well as obtaining the 
dynamical evolution through a BSE instead of a SDE [see, for example, Eq. (14. 5p . Eq. (I8.25P 
and Eq. <K2M . 



43 



A central result of the present work is the formal equivalence between the massless bound- 
state formalism and the corresponding approach based on the direct study of the SDE of 
the gluon propagator 2l|. In particular, the powerful relations provided by the STIs of the 
theory, such as Eq. (I8.18p . allowed us to demonstrate the exact coincidence of the integral 
equations governing the momentum evolution of the gluon mass in both formalisms. 

Note that this formal equivalence, in addition to serving as a clear indication of an under- 
lying consistency, opens up the possibility of extracting useful, albeit indirect, information 
on the structure of certain quantities appearing in one formalism, from results obtained 
within the other. It seems, for instance, that a nonperturbative approximation for the quan- 
tity Y (see the SDE of Eq. (18.2 ip ) is easier to obtain than a corresponding approximation 
for the four-gluon kernel /C, a central ingredient of the BSE in Eq. (I8.24p . Indeed, Y involves 



the fully dressed three-gluon vertex, which may be partial 



gauge-technique Ansatz, from the STIs that it satisfies 56]; to be sure, lattice simulation 



y reconstructed, by means of a 



may also be valuable in this effort [69j. Instead, nonperturbative approximations for /C are 
much more difficult to obtain with the present technology. Therefore, one might be able to 
use a finer solution for m 2 (q 2 ), obtained from the SDE of Eq. (I8.2ip by going beyond the 
one-loop approximation for Y, to infer the behavior of the kernel /C, at least for some special 
kinematic configurations. This information, in turn, may be helpful in a variety of unrelated 
studies that involve this particular kernel. 

Even though the derivation of various of the results presented here relies on the use 
of the Landau gauge, it would seem that the general features of the massless bound-state 
formalism remain valid for any value of the gauge-fixing parameter, within the general class 
of linear covariant gauges. It would be therefore most interesting to derive the corresponding 
dynamical equations for different gauges, and explore whether the gluon mass generation 
mechanism persists away from the Landau gauge. In fact, such studies may be com plem ented 
by parallel lattice simulations of the corresponding gluon and ghost propagators |70l. \71\. 

In a similar vein, lattice simulations carried out in the BFM 72|, [73] may prove instru- 
mental for verifying explicitly some of the formal relations employed throughout this work. 
For example, the BQI of Eq. (15. ip may be directly probed, for a wide range of momenta, if 
A(g 2 ) is simulated in the background Landau gauge; for a prediction, see 64|. In fact, even 
the knowledge of just the corresponding saturation point, A(0), would enable one to check 
the validity of Eq. (ETJ) at q 2 = 0. 
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